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FOREWORD 


This  work  was  performed  under  US  Army  Natick  Laboratories  Contract 
No.  DAAG17 -70-0-0127  during  the  period  of  1  Apr  70  to  31  Mar  71 »  The 
Project  No.  was  1F162203D195>  entitled  "Exploratory  Development  of 
Airdrop  Systems",  and  the  Task  was  No.  13  entitled  "Impact  Phenomena" „ 
Messrs.  Edward  J .  Giebutowski  and  Marshall  S .  Gustin  of  the  Airdrop 
Engineering  Laboratory  served  as  the  Project  Officers. 

The  effort  is  part  of  a  continuing  investigation  directed  towara 
obtaining  a  better  understanding  of  the  failure  mechanism  of  energy 
dissipater  materials,  and  the  response  of  airdroppable  supplies  and 
equipment  to  airdrop  impact  phenomena;  and  toward  obtaining  improved 
airdrop  energy  disipater  materials  and  techniques. 

This  report  is  concerned  primarily  with  an  analytical  study  of  the 
configuration  change  of  an  individual  cell  of  paper  honeycomb  energy 
dissipater  material  as  it  is  crushed,  in  order  to  determine  the  relative 
effect  of  changes  in  various  cell  design  characteristics  on  level  of 
strength  and  uniformity  of  performance. 
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is.  abstract  jKrcflgm&ra cajoy  exact  equations  or  tne  aenections  required  lor  a  circular 
cylinder  to  buckle  into  a  developable  polyhedral  shape  are  derived.  From  this 
equation  it  i3  seen  that  the  tangential  displacement  can  be  readily  related  to 
the  radial  displacement.  The  exact  formulas  for  coefficients  of  a  Fourier  series 
representing  the  radial  displacements  are  derived  for  a  buckled  circular  cylinder 
having  ary  longitudinal  and  circumferential  mode  numbers.  This  is  in  contrast  to 
earlier  work  in  which  approximate  Fourier  coefficients  were  derived  for  large  values 
of  circumferential  mode  number,  (n£lO). 

The  analytical  expressions  are  derived  to  study  buckling  of  individual 
hexagonal  cells  in  paper  honeycomb.  The  honeycomb  cells  are  originally  in  the 
developable  shape  represented  by  the  circumferential  mode  number  n*6  and  longitlidinal 
mode  number  m"0*  With  the  guidance  of  the  analytical  expressions  described  above, 
it  has  been  observed  that  individual  honeycomb  cells  constrained  by  surrounding 
honeycomb  cells  in  an  impacted  honeycomb  pad,  buckle  in  a  modified  form  of  the  n*=3 
deflection  shape.  This  modified  deflection  shape  is  determined  as  a  function  of  the 
width  of  the  glue  line  joining  the  cells  which  remain  intact  after  crushing,  that 
is  the  effective  glue  line  width. 

The  analytical  expressions  of  the  Fourier  coefficients  of  the  radial 
displacements  of  the  modified  n*3  deflection  shape  are  derived.  The  net  radial 
displacements  of  the  impacted  buckled  hexagonal  shell  ax*  obtained  by  subtracting 
the  Fourier  series  representing  the  original  hexagonal  cylindrical  shape  corresponding 
to  the  n=6,  m*0  mode  from  the  Fourier  series  representing  the  modified  n=3  buckling 
shape.  This  net  deflection  is  inserted  in  the  strain  energy  equations  to  obtain 
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the  total  elastic  strain  energy,  bending  and  stretching  strain  energies, 
and  buckling  stress.  These  quantities  are  obtained  for  hexagonal  shells 
with  variations  in  various  geometrical  parameters,  such  as  radius/ thickness, 
radius/length,  effective  glue  line  width  and  the  mechanical  properties; 
Poisson's  ratio  and  Young's  modulus. 
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ABSTRACT 


To  obtain  some  insight  into  the  failure  mechanism 
of  paper  honeycomb  and  to  identify  the  significant  para¬ 
meters  an  analysis  has  been  made  of  the  buckling  of  a 
thin  walled,  axially  loaded  cylinder  with  a  hexagonal 
cross  section.  The  analysis  begins  with  a  derivation  of 
mathematically  exact  expressions  for  the  radial  deflec¬ 
tions  required  to  change  a  circularly  cylindrical  shell 
into  a  developable  polyhedral  shape.  This  expression  is 
then  put  into  the  form  of  a  Fourier  series  which  repre¬ 
sents  the  radial  displacements  of  a  buckled  circular 
cylindrical  shell  having  any  longitudinal  and  circumfer¬ 
ential  mode  numbers.  Hexagonal  honeycomb  cells  have  the 
developable  shape  given  by  that  Fourier  series  for  the 
circumferential  mode  n  =  6  and  the  longitudinal  mode 
m  =  0.  From  this  analytical  expression  it  is  observed 
that  individual  honeycomb  cells  in  an  impact  loaded 
honeycomb  pad  buckle  in  a  modified  form  of  the  n  =  3 
deflection  shape.  This  modified  shape  is  determined  as 
a  function  of  the  width  of  the  glue  line  joining  the 
cells.  This  is  referred  to  as  the  effective  glue  line 
width. 

The  net  radial  displacements  of  the  impact  buckled 
hexagonal  cell  are  obtained  by  subtracting  the  Fourier 
series  representing  the  original  hexagonal  cylindrical 
shape  corresponding  to  the  n  =  6 ,  m  =  0  mode  from  the 
Fourier  series  for  the  modified  n  =  3  buckling  shape. 
This  net  deflection  is  inserted  in  the  strain  energy 
expressions  for  the  total  elastic  strain  energy,  inclu¬ 
ding  both  bending  and  stretching  strains.  Using  these 
energy  expressions  the  buckling  stress  is  computed. 

Buckling  stresses  are  determined  for  hexagonal 
cells  in  which  variations  in  the  geometrical  parameters, 
radius/thickness,  radius/length,  and  effective  glue 
line  width  are  introduced.  In  addition  the  mechanical 
properties,  Poisson's  ratio  and  Young's  modulus  are 
varied. 

The  effective  glue  line  width  is  shown  to  be  the 
most  significant  parameter. 
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INTRODUCTION 


1 .  Background 

Paper  honeycomb,  shown  in  Fig.  1,  has  bee/i  used  for 
the  past  15  years  as  an  energy  dissipation  device  for 
cushioning  supplies  and  equipment  airdropped  by  parachute. 

Its  outstanding  advantages  are  low  cost,  lightweight,  and 
a  crushing  strength  that  is  essentially  independent  of 
variations  in  parameters  such  as  impact  velocity  and  mois¬ 
ture  content. 

The  configuration  of  the  honeycomb  pads  presently 
used  for  impact  cushioning  was  chosen  on  an  empiricaJ  basis, 
because,  unfortunately,  the  critical  parameters  in  the 
energy  dissipating  process  have  not  yet  been  identified. 

It  is  known  that  the  crushing  stress  increases  with  the 
honeycomb  density,  but  how  the  cell  size,  paper  weight, 
and  glue  line  width  individually  affect  the  crushing  stress 
has  not  been  determined.  One  of  the  particularly  puzzlinc? 
aspects  of  the  behavior  of  paper  honeycomb  has  been  the 
variation  observed  in  the  crushing  strength  of  apparently 
identical  honeycomb  pad  specimens. 

The  purpose  of  this  investigation  is  to  study  analy¬ 
tically,  the  effects  of  cell  size,  wall  thickness,  and 
glue  line  width  on  the  crushing  load  of  a  honeycomb  pad. 

The  primary  objective  of  the  study  has  been  the  development 
of  a  rational  guide  for  the  design  and  construction  of 
honeycomb  pads  for  various  specific  applications.  It  was 
also  intended  that  the  analysis  would  lead  to  energy  dissi¬ 
pation  configurations  superior  to  those  of  honeycomb  presently 
in  use. 

Thi  approach  taken  in  this  investigation  is  primarily 
analytical,  guided  by  limited  experimentation.  An  analy¬ 
tical  rather  than  experimental  approach  was  taken  because 
of  the  difficulty  in  obtaining  consistent  test  results 
even  with  specimens  produced  by  a  rigidly  controlled  fabri¬ 
cation  process  in  which  no  parameter  variations  are  included. 

In  order  to  ach: eve  this,  finite  deflection  buckling 
of  a  single  cell  is  analyzed  as  though  the  stresses  stay 
within  the  elastic  range  during  the  entire  buckling  process. 
Actually  it  is  not  unreasonable  to  assume  elastic  buckling. 
Plastic  deformation  has  to  occur  for  significant  energy 
dissipation,  but  it  occurs  for  the  most  part  after  the 
initial  buckling  has  occurred. 

The  buckle  pattern  of  cells  in  honeycomb  pads  crushed 
statically  or  dynamically  resembles  the  buckle  pattern  of 
circular  cylindrical  shells  subjected  to  axial  impact  as 
shown  in  Fig.  2b.  the  original  approach  to  this  problem 
it  was  intended  tha  •  critical  parameters  would  be  determined 
by  representing  the  honeycomb  cells,  which  initially  have 
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Pig.  2  Ideal  Inextensional  and  Actual  Collapse  Patterns . (Prom  Cuppa) 


rounded  hexagonal  cross  section  shapes,  as  circular 
cylindrical  shells  subjected  to  impact  loading.  Tfie  disr- 
covery  of  a  Fourier  series'  expansion  for  the  developable 
surfaces  made  it  possible  to  investigate  the  behavior  of 
a  hexagonal  cell.  s  ,  i 

2.  Previous  Work  > 

Much  of  the  most  significant  analytical  and  experi¬ 
mental  work  concerning  longitudinal  impact  'loading  of  thi’n  > 
circular  cylindrical  and  conical  shells  has  been  done  by 
Coppa.  He  found  that  the. buckled  sljiape  of  impacted 

cylindrical  shells  is  similar  .to  the  developable  polyhedral 
shape  qualitatively  described  by  Yoshimura4  and ; shown  in 
Fig.  3. 

.1  •  ‘  i 

a.  Shell  Impact 

Coppa  suggested  that  buckling  of  cylindrical  and  coni¬ 
cal  shells  subjected  to  longitudinal  impact  is  a1  phenomenon 
related  to  the  propagation  of  stress  v^aves  generated  at :  the 
impacted  end  of  the  shell.  •  .  .  ! 

The  induced  stress ’wave  propagates  through  ;the  length 
of  the  shell  until  it  exceeds  the  critical  buckling  stress 
of  a  section  of  the  shell,  length..  If  the  initial  stress  . 
produced  at  the. impacted  end  is  not  large  enough  to  exceed 
the  critical  buckling  stress,  the  initial =stress  amplitude 
is  successively  ipcreased  by  reflection  at  the  ends  of  the 
shell  during  stress  wave  propagation  until  a  critical  local 
buckling  stress  is  reached  and  local  buckling  commences.  ! 
This. stress  increase  phenomenon  is  produced  during  stress 
wave  propagation  by  reflection  from  £he  higher  acoustic  > 
impedance  boundaries  at  the  shell  ends.  The  acpustic 
impedance  is  defined  as  the  'product  of  the  square  root  '  1 

of  the  elastic  modulus  and  the  niass  density  ef  a  material.' 
Coppa* s  hypothesis  is  corroborated  to  some  extent  by  the 
experimental  observation  that  in  his  cylinders  and  in 
University  of  Texas  honeycomb  test  specimens,  the  buckles 
tend  to  occur  first  at  the  top’  and  bottom  of  ‘the  shell  or 
cell  walls.1  Honeycomb  cells  appear  to  buckle  at  'midlength 
only  after  very  extensive  crushing.  This  characteristic 
may  be  due  to  the  buckle  inhibiting  effect ?  of  entrapped  air. 

b.  Developable  -Polyhedral  Surface  ' 

Yoshimura  pointed  out  that  it  is  possible  for  a  cylin¬ 
drical  surface  to  buckle  into  a  developable  polyhedral 
shape  as  shown  in  Fig.  3b.  This  developable  polyhedral 
shape  is  composed  of  a  number  of  plane  surface  panels  which 
have  no  extensional  strain.  There  is  infinite  curvature 

at  the  shard  fold  edges.  Yoshimur.a  suggested  that  the  exact 
deflections  for  this  developable  surface  might  te'  represented 
by  a  Fourier  series  with  an  infinite  number  of  terms. 

Furthermore,  Yoshimura  suggested  that  the  deflections 
of  buckled,  finite  thickness,  cylinders  with  the  small  radii 
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of  curvatures  that  occur  at  the  rounded  fold  edges  as  shown 
in  Pig.  2b,  could  be  approximated  by  a  truncated  form  of 
the  infinite  Fourier  series.  Yoshimura  used  only  a  few 
terms  in  his  Fourier  series  representation  of  the  displace¬ 
ments  required  for  the  developable  surface  and  described 
the  surface  more  in  qualitative  than  in  quantitive  terms. 

■  Hoff,  Madsen  and  Mayer5  derived  an  approximate  formula 
for  the  coefficients  of  a  Fourier  series  representing  the 
deflections  of  a  cylindrical  surface  deformed  into  a  poly¬ 
hedral  surface.  This  formula  is  apparently  valid  for 
values  of  the  circumferential  mode  number  n  equal  to  or 
greater  than  10.  Yoshimura’ s  hypothesis  of  the  developable 
shape  was  verified  by  comparing  the  values  of  the  coeffi¬ 
cients  of  a  Fourier  series  for  describing  the  buckled  shape 
of  very  thin  axially  loaded  cylinders.  It  wjs  found  that 
the  larger  the  number  of  terms  used  in  the  series,  the 
closer  the  shape  approached  the  developable  shape  qualita¬ 
tively  described  by  Yoshimura. 

3.  Present  Investigation 

A  few  static  and  impact  tests  were  performed  on  shells 
in  order  to  observe  the  buckling  patterns  and  crushing 
strengths  under  both  types  of  loading.  Observation  of  the 
buckle  patterns  was  helpful  in  the  derivation  of  formulas 
for  the  displacements  required  to  produce  developable 
surfaces. 

The  equations  for  the  exact  radial  displacements 
required  for  a  circular  cylindrical  surface  to  deform  into 
a  developable  polyhedral  surface  are  derived  from  geometri¬ 
cal  considerations.  A  formula  for  the  coefficients  of  a 
Fourier  series  expansion  for  the  exact  radial  displacements 
valid  for  all  values  of  n  is  then  derived.  A  simple  relation 
which  exists  between  the  radial  and  tangential  displacements 
eliminates  the  necessity  of  determining  Fourier  coefficients 
for  the  tangential  displacement.  '  It  Is  shown  that  the  ini¬ 
tial  hexagonal  cross  section  can  be  expressed  in  a  Fourier 
series  representation,  which  provides  a  better  model  than 
the  circular  "croSS* ^detlon  shell  ‘Initially  assumed.  The 
buckled  shape  of  honeycomb  cells  is  found  to  be  a  modified 
form  of  a  developable  surface  into  which  a  circular  cylin¬ 
der  can  deform.  This  modified  shape  is  a  function  of  the 
effective- glue  line  width  between  adjacent  honeycomb  .cells 
in  the  honeycomb  pad. '  Details  of  these  derivations  are 
given  in  the  Appendix. 

Once  analytical  expressions  for  the  deflections  in  the 
buckled  shape  are  available,  an  energy  method  can  be  used 
for  computing  buckling  loads. 
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PROCEDURE 


A  single  honeycomb  cell  will  be  considered  because  it 
is  simpler  to  study  analytically  the  response  of  single 
cells  than  it  is  to  study  clusters  of  cells.-  Furthermore, 
individual  cells  may  be  responding  differently  depending 
on  their  physical  location  in  the  honeycomb  composite. 

In  experimental  investigations  it  is  simpler  to  study  the 
behavior  of  groups  or  clusters  of  cells.  The  buckled  shape 
of  a  honeycomb  cell  resembles  the  developable  polyhedral 
surface  which  has  slope  discontinuities  or  zero  radii  of 
curvature  as  shown  in  Fig.  3b.  The  actual  buckled  cylinder 
with  a  finite  thickness  has  folds  which  are  more  rounded. 

The  roundness  of  the  corners  at  the  folds  increases  as  the 
cylinder  thickness  increases. 

Theoretically,  a  function  with  slope  discontinuities 
or  infinite  curvature  can  be  exactly  represented  by  a 
Fourier  series  \*ith  an  infinite  number  of  terms  (constrained 
only  by  the  Dirichlet  conditions) . 6  if  this  Fourier  series 
is  truncated  it  will  represent  a  surface  with  a  finite 
curvature  at  the  position  of  the  original  slope  discontin¬ 
uities.  The  curvature  at  the  folds  increases  with  the  num¬ 
ber  of  terms  retained  in  the  series. 

If  the  exact  displacement  equations  could  be  repre¬ 
sented  identically  by  a  Fourier  series  with  an  infinite 
number  of  terms,  there  would  be  zero  extensional  strain 
energy  but  infinite  bending  moment  for  elastic  bending 
at  the  fold  edges  for  the  model  shown  in  Fig.  3b.  This 
Fourier  series  when  truncated  to  a  finite  number  of  terms, 
should  provide  a  good  representation  of  the  deflection 
shape  of  a  buckled  finite  thickness  cylinder.  An  increasing 
number  of  terms  would  be  required  as  the  cylinder  wall 
becomes  thinner.  This  is  because  the  thinner  the  buckled 
cylinder,  the  closer  it  approaches  the  developable  shape 
assumed  by  a  buckled  infinitely  thin  circular  cylinder. 
Increasing  the  number  of  Fourier  series  terms  used  would 
also  decrease  the  extensional  strain  energy  and  increase 
the  bending  strain  energy.  For  an  infinite  number  of 
terms  the  extensional  strain  energy  vanishes  and  the 
bending  strain  energy  approaches  infinity. 

One  way  to  obtain  the  buckling  load  in  a  static 
analysis  is  to  set  the  external  work  done  on  the  cylinder 
equal  to  the  total  strain  energy  in  the  cylinder.  This 
will  in  general  not  be  the  minimum  buckling  load.  To  obtain 
the  minimum  the  number  of  terms  in  the  Fourier  series  for 
the  deflections  is  varied  until  a  sufficient  number  of 
buckling  loads  has  been  computed  to  establish  the  minimum. 
The  buckled  shape  is  determined  by  the  number  of  terms 
included  in  the  series. 
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1.  Developable  Surface  Deflection  Shape 

Since  the  buckled  shape  of  a  cylinder  has  symmetry  and 
periodicity  in  the  circumferential  coordinate  0  and  the 
axial  coordinate  x,  as  shown  in  Fig.  4,  only  the  deflections 
of  one  panel  are  necessary  to  obtain  the  coefficients  for 
the  Fourier  series  representation  of  the  displacements  when 
suitable  coordinate  transformations  are  made.  The  deriva¬ 
tions  of  the  exact  radial  and  tangential  displacement  expres¬ 
sions  are  given  in  Section  1  of  The  Appendix.  Eqs. (l.llj 
and  (1.12),  the  expressions  for  the  radial  and  tangential 
displacements  necessary  to  transform  a  Circular  cylindrical 
shell  into  a  developable  polyhedral  surface  show  that  there 
is  a  simple  relation  between  the  tangential  displacement 
and  the  first  derivative  of  the  radial  displacement.  This 
allows  the  tangential  displacement  to  be  expressed  approxi¬ 
mately  in  terms  of  the  Fourier  series  expansion  of  the 
radial  displacement.  Then  there  is  no  necessity  for  deter¬ 
mining  the  Fourier  coefficients  for  the  tangential  displace¬ 
ment. 

The  radial  displacement  functions  are  expanded  in  a 
double  Fourier  series  in  0  and  x  in  the  following  manner. 
First,  the  two  derived  functions  of  the  exact  radial  dis¬ 
placement  w,  with  limits  as  a  function  of  x,  are  represented 
by  a  single  Fourier  series  expansion  in  the  circumferential 
coordinate  0 .  This  Fourier  series  has  coefficients  a  which 
are  a  function  of  x.  The  formula  for  the  a  coefficients 
is  determined  by  integrating  the  derived  exact  radial  dis¬ 
placements  over  two  regions  of  0  which  have  limits  that 
depend  on  x. 
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wKe)  *=  ^.-La-jet) cos jn& 

The  deflected  shape  obtained  has  slope  discontinuities 
in  the  axial  direction,  but  with  a  finite  number  of  terms 
in  this  expansion  there  is  no  slope  discontinuity  m  the 
circumferential  direction.  These  a  (x)  Fourier  coefficients 
are  themselves  now  represented  by  a^ Fourier  series  in  terms 
of  x,  such  that 
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The  coefficients  a,,  in  this  Fourier  series  are 
derived  as  shown  in  theJ Appendix.  When  an  infinite  number 
of  terms  are  taken  in  this  Fourier  series  it  yields  exact 
displacement  values  and  slope  discontinuities  as  shown  in 
Fig.  3b.  With  a  finite  number  of  terms  in  the  series  there 
are  no  slope  discontinuities  in  either  the  circumferential 
or  axial  directions. 

A  computer  program  to  compute  the  Fourier  coefficients 
a,(x)  and  a.,  was  written  for  the  CDC  6600  computer.  The 
computer  was'' used  to  determine  the  coefficients,  plot  the 
surface  deformations  and  compare  the  Fourier  series  repre¬ 
sentation  of  the  radial  deflection  with  the  exact  values 
of  the  radial  deflection. 

Hoff,  Madsen  and  Mayer5  have  published  formulas  for 
computing  approximate  values  of  the  Fourier  coefficients 
for  large  circumferential  mode  numbers  (n  ^  10).  A  compu¬ 
ter  program  was  written  in  order  to  obtain  the  values  of 
these  approximate  coefficients.  These  approximate  coeffi¬ 
cients  are  compared  in  RESULTS,  with  the  more  exact  values 
of  aij  referred  to  above. 

2.  Honeycomb  Deflection  Shape 

After  the  coefficients  of  the  Fourier  series  represen¬ 
tation  of  the  radial  displacement  for  the  developable 
shape  are  obtained,  it  is  possible  to  determine  the  radial 
displacements  required  to  produce  the  original  hexagonal 
cylindrical  shell  shape  of  individual  cells  of  the  honeycomb 
pad.  The  hexagonal  cell  is  merely  the  n=6,  m=0  mode  of  the 
circular  cylindrical  shell  as  shown  in  Fig.  5.  This  shape 
represents  the  nearly  perfect  hexagonal  shape  of  precision 
made  honeycomb.  In  mass  produced  manufactured  honeycomb, 
the  cell  tends  toward  a  rounded  hexagonal  shape,  which  can 
be  more  closely  represented  by  a  smaller  number  of  terms 
in  the  Fourier  series  for  the  n=6,  m=0  mode. 

During  this  work,  to  establish  the  nature  of  the 
buckling  pattern  in  crushed  honeycomb,  an  impact  test  was 
performed  in  which  the  honeycomb  pad  was  cut  into  four 
pieces  prior  to  testing.  The  cuts  were  made  along  the  cell 
edges  and  then  the  four  cut  pieces  were  butted  together  for 
the  test.  This  approach  was  used  to  avoid  the  distortion 
in  the  buckle  pattern  which  cutting  after  testing  usually 
produces.  Three  of  the  pieces  illustrating  the  test  config¬ 
uration  are  shown  in  Fig.  1.  The  cell  buckle  pattern  is 
shown  in  Fig.  6. 

From  this  photograph  it  may  be  observed  that  the  buckling 
pattern  of  the  cells  in  the  crushed  honeycomb  is  such  that 
the  radial  displacement  from  the  sides  which  formed  the 
original  hexagonal  cross  section  shape  alternates  from  inward 
to  outward  as  one  progresses  around  the  periphery  of  the  cell 
wall  and  the  maximum  displacements  appear  to  be  equal. 
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As  mentioned  previously,  the  initial  honeycomb  hexa¬ 
gonal  shape  is  equivalent  to  the  n «6,  m-0  mode  of  a  circular 
cylinder.  Cross  sections  of  a  buckled  cylindrical  shell  in 
the  n=3  mode  have  three  sides  which  form  an  equilateral  tri¬ 
angle,  at  the  axial  nodes.  However,  there  are  six  sides  at 
any  cross  section  between  the  axial  nodes.  Maximum  displace¬ 
ments  occur  at  the  circumferential  midpoints  of  the  original 
sides  of  the  hexagon  as  shown  in  Pig.  7.  Prom  the  exact 
equations  for  the  radial  displacements  it  is  found  that 
the  difference  between  the  maximum  inward  deflection  at 
these  points  for  the  n=3  mode  and  the  n-6,  m*0  mode,  is 
equal  to  the  maximum  outward  deflection.  Thus  the  alter¬ 
nating  inward  and  outward  displacements  observed  in  Fig.  6 
can  be  regarded  as  the  displacements  of  a  circular  cylin¬ 
drical  surface  in  a  modified  n=3  buckling  mode  minus  the 
n-6,  m=0  mode  which  represents  the  original  hexagonal  cylin¬ 
drical  shape  of  the  cell  shown  in  Fig.  5. 

The  shape  of  the  modified  n=3  buckling  mode  of  the 
honeycomb  cells  is  found  to  be  a  function  of  the  width  of 
the  glue  line  bonding  adjacent  cell  walls  which  remains 
intact  during  buckling,  that  is,  the  ’'effective  glue  line 
width."  The  cross  section  of  a  buckled  circular  cylindrical 
shell  at  the  midpoint  (1/2  in  Pig.  7)  between  the  axial 
nodes  of  the  n=3  buckling  mode  is  a  regular  hexagon.  If 
there  is  no  extension  of  surfa.is  this  section  for  the  n=3 
mode  must  be  the  same  as  the  cross  section  for  the  n*6,  m=0 
mode.  Thus  there  is  no  change  from  the  original  hexagonal 
shape.  The  axial  location  of  this  section  corresponds  to 
the  axial  nodes  of  a  buckled  honeycomb  pad  cell.  The  ends 
of  the  honeycomb  cells  where  displacements  are  constrained 
to  the  hexagonal  cross  section  shape  by  the  facing  paper 
will  always  be  axial  nodes. 

It  is  very  difficult  to  tell  by  examination  of  buckled 
honeycomb  exactly  what  the  geometry  of  the  buckled  cell  is 
Clearly  there  are  alternating  inward  and  outward  radial 
deflections.  However,  the  hexagonal  cells  are  restrained 
at  4  points,  i.e.,  there  are  ^  glue  lines,  not  three  as 
shown  in  Fig.  8.  Two  of  the  four  glue  lines  act  as  point 
restraints  rather  than  lines.  Thus  it  seems  very  certain 
that  the  buckling  pattern  cannot  be  identical  to  the  one 
for  the  nrs3,  or  the  modified  n=3  modes.  Nevertheless, 
these  shapes  which  can  be  represented  in  a  reasonably 
convenient  mathematical  form  can  be  used  to  approximate 
the  buckling  behavior  of  honeycomb  cells  in  a  pad,  and  to 
obtain  an  indication  of  the  effects  of  the  restraint  offered 
by  the  glue  lines. 

To  continue  the  analysis  of  the  glue  line  effects  it  is 
necessary  now  to  express  the  radial  deflections  for  the  modi¬ 
fied  n=3  mode  in  an  analytical  form  that  can  later  be  used 
in  strain  energy  calculations.  The  necessary  equation  can 
be  obtained  by  taking  the  equation  given  previously  for 
,  i.e.. 
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Fig.  8  Profile  of  Circular  Cylinder  Buckled  in  n=6,  m=0;  n=3  and  Modified  n=3  Modes 


and  computing  the  coefficients  a. -  for  the  deflections  shown 
in  Fig.  8c  for  the  modified  n=3  mode.  These  deflections 
are  obtained  by  applying  purely  geometrical  considerations 
and  the  requirement  of  no  extensions  in  a  developable  surface. 

When  n  is  specified,  the  number  of  folds  in  the  circum¬ 
ferential  directions  is  specified.  This  does  not  mean  that 
the  paper  honeycomb  will  necessarily  buckle  in  this  pattern. 
Consequently  if  n  is  specified  all  conclusions  regarding 
buckling  loads  must  be  predicated  on  the  assumption  that 
the  actual  buckling  pattern  will  agree  with  the  assumed 
pattern.  The  number  of  longitudinal  folds  (buckles)  is 
not  specified.  The  more  folds  or  buckles  there  are  the 
higher  the  buckling  load  is  likely  to  be. 

It  will  now  be  shown  that  if  a  cell  buckles  into  the 
modified  n=3  mode  the  effect  of  the  glue  line  will  be  to 
decrease  the  buckle  length,  and  increase  the  number  of  longi¬ 
tudinal  folds,  and  hence  the  buckling  load. 

Consider  Fig.  8a.  This  is  one  of  the  polyhedral 
panels  of  the  buckled  cell  surface  viewed  normal  to  the 
surface.  The  solid  lines  indicate  a  panel  for  modified 
n=3  buckling,  and  the  dotted  lines  indicate  the  shape 
for  the  plain  n=3  mode.  For  this  mode  the  panel  is  an 
equilateral  triangle  with  sides  277a/3.  (The  inextensibility 
conditions  requires  the  total  length  of  three  sides  to  be 
the  same  as  the  circumference  of  the  cylinder,  diameter  2a, 
from  which  the  developable  polyhedral  surface  is  formed.) 

For  the  modified  n=3  mode  the  triangular  panels  are  changed 
to  trapezoids.  If  the  glue  line  width  is  2b,  as  seen  in 
Fig.  3b,  the  sides  of  the  trapezoid  have  a  length 
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as  shown  in  Fig.  8a.'  The  line  AB  in  Fig.  8c  is  the  line 

AB  shown  in  Fig.  8a.  The  projection  of  AB  on  the  horizontal 

in  Fig.  8c  is  A  and  where  L  is  the  initial  length  of 

the  cylinder. .  The  horizontal  projections  of  AB  and  A'B'  in 
Fig.  8c  are  in  the  same  proportion  as  the  sides  of  the  tri¬ 
angle  and  the  trapezoid  in  Fig.  8a.  Hence 


or  A  2 
fied  n=3  mode. 


(iiro  _  = 


TT0-}**'  =  the  buckle  length  of  the  modi- 

The  projection  of  AA'  cn  the  horizontal  is 
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the  same  fractional  Dart  of  as  2b  is  of  2Tra  ,  that  is, 
as  shown  in  Fig.  8c.  2 

To  represent  the  modified  n=3  mode  the  double  Fourier 
series  for  the  radial  deflection  can  be  obtained  from  the 
equation  for  the  n=3  mode  by  simply  modifying  the  argument 
of  the  cos  imiTx  term  to  take  into  account  the  change  in 

the  Deriodic  interval  which  the  glue  line  causes.  Thus 
for  the  modified  n=3  mode 

oo  oo 

w3(x,q)=  2  2  *'J  cos FT^i)p  cos  3Je 
*•!  j wi'  * 

where  X,  =  X  -  Jf 

The  coordinate  transformation  is  introduced  to  out 
the  origin  of  the  coordinate  system  at  point  A  ' ,  the  fold 
in  the  buckling  pattern. 

The  complete  derivation  of  the  exDressior.s  for  the 
a^  coefficients  is  given  in  section  4  of  the  Aopendix. 

3.  Buckling  Stress 

The  Fourier  series  exoansion  for  the  net  deflection 
of  the  buckled  honeycomb  cell  is  inserted  in  the  bending 
and  stretching  strain  energy  equations  as  shown  in  Section 
5  of  the  Appendix.  To  obtain  a  buckling  stress  the  sum  of 
the  bending  and  stretching  strain  energies  is  equated  to 
the  work  done  by  the  aDnlied  axial  force  ud  to  the  time  the 
buckling  pattern  is  formed.  These  quantities  are  each 
functions  of  the  deflection  shape.  The  deflection  shape  is 
also  a  function  of  the  indices  i  and  j  in  the  double  Four¬ 
ier  series.  The  higher  the  values  of  i  and  J  the  closer  the 
developable  shape  is  aoproached.  The  indices  are  increased 
until  a  minimum  buckling  stress  is  obtained.  This  value 
represents  the  stress  required  to  buckle  the  shell  into  a 
Darticular  mode.  (In  this  case  the  modified  n=3  mode.) 

The  non-dimensional  buckling  stress,  o/E  is  shown 
as  a  function  of  geometrical  parameters  t/c,  1/ c,  the  glue 
line  width,  2b,  and  Poisson's  ratio  in  the  RESULTS. 
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RESULTS 


1.  Developable  Polyhedral  Surface  Equations — Exact  Displacement 
The  shape  of  a  circular  cylinder  which  has  buckled  into 
a  polyhedral  developable  surface  is  symmetric  and  periodic  i 
in  the  circumferential  coordinate  and  the  axia]  coordinate  x. 
Thus  only  the  eauations  of  the  deflections • of  a  single  .nanel 
of  the  buckled  surface  are  Initially  required  to  obtain  the 
Fourier  series  expansion  of  .the  displacement  of  the  entire 
buckled  surface  when  suitable  coordinate  transformations  are 
made.  . 

The  radial  and  circumferential  'displacements,  v:  and  v 
respectively,  of  a  single  panel  have  each  been  determined  :  ; 

purely  geometrical  analysis  as  two  separate  functions  c'# 
the  circumferential  coordinate  0  and  axial  coordinate  x. 

Each  function  is  valid  in  a  region  whose  .c  Irciunfer'ont  ’  a i 
coordinate  limits  are  a  function  of  the  ax.iai  coord’ rate-  x.  1 
These  functions  derived  in  flection  1  of  the- Aru-or/':’ x  ar"-:. 

1  (  * 

W  (x,&)  -  [a.  +  c2(x)J  cosG  '+  slQ  sin  9  -a.  , 

.  os  e*  [/-£]$ 

V  ( X,  6)  -  [a.  +  cgM] sin  O  -  a.  Geos  O 


*  i 

W  (x,9)  =  [a.  +  c,{x)]co$f2-e)+(%-6)sis/n(%-G)-a. 

,  tO-th***  , 

V  (x,0)  ={*  *  C,  (x)]'sin  (*-e)+(3-G)AC0s(%-e) 


where 


o  -  L 

c*  (x)  —  (  -  )x  ~  A  =  Radi*/  Pispftcemenf amplitude  * 

x  at  xf  9  =0 
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and  Q  and  A  are  the  amplitudes  of  the  maximum  radial 
displacement  at  the  longitudinal  nodes  in  the  outward 
and  inward  directions  respectively.  The  values  of  the 
amplitudes  given  by  Coppa^  are: 

<?  -  sari* 


From  Eqs.  (1)  and  (2)  it  may  be  seen  that  there  is 
a  simple  relationship  between  the  tangential  and  radial 
displacements.  The  tangential  displacements  can  be 
expressed  in  the  following  manner: 


i r(x,0) 
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a.  Fourier  Coefficients  of  Double  Fourier  Series 

The  derived  radial  displacement  functions  w  given  in 
Eq.  (1)  are  expanded  into  a  double  Fourier  series  in  •? 
and  x  as  follows.  Complete  details  of  the  expansion  are 
given  in  section  2  and  3  of  the  Apoendix. 

First  the  two  functions  of  radial  displacement  w 
given  in  Eq.  (1)  are  represented  by  a  single  Fourier  series 
expanded  in  9 .  The  Fourier  coefficients  a .  are  a  function 
of  x.  '  3 

as 

*/(x)9)  =  £  eij(x)COSjn0 


The  deflected  shape  obtained  has  slope  discontinuities  in 
the  axial  direction/  but  with  a  finite  number  of  terms  in 
this  expansion  there  is  no  slope  discontinuity  in  the  cir¬ 
cumferential  direction. 

The  a . (x)  coefficients  are  determined  by  integrating 
the  two  functions  of  w  from  Eqs.  (1)  ans  (2)  over  the  two" 
applicable  regions  of  e(x)  as  presented  in  Section  2  of 
the  Appendix. The  equation  derived  is: 
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where 


e  =  /  *jn 
F  =  l-jn 

a  —  m  Try. 

~n7T 

These  a . (x)  Fourier  coefficients  are  expanded  in  a 

Fourier  series  in  terms  of  x  in  Section  3  of  the  Appendix  such 

that  „ 

oo 
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or 
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atJ  cos 
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w  (x,  6)  =  £  a,jj  cos  -  cos jnO 
J'°  i-0 


An  infinite  number  of  terms  in  this  series  gives  the 
deflection  of  a  developable  surface  with  slope  discontinuities 
at  the  surface  folds.  However  a  finite  number  of  terms  repre¬ 
sents  folds  with  finite  radii  of  curvature. 

The  formula  for  the  Fourier  coefficient  a. .  is  obtained 
by  integrating  the  Fourier  coefficient  a .  (x)  in-3  the  region 
of  one  axial  half  wave  length  by  taking  Into  consideration 
the  symmetry  of  the  deflections  in  the  axial  direction. 

This  derivation  is  given  in  detail  in  Section  3  of  the  Anoendix, 
The  equation  is 
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2.  Computer  Check 

Numerical  values  of  a,.  have  been  computed  and  used 
to  plot  the  profile  of  the  deflections  of  a  circular  cylinder 
buckled  in  the  n=3  mode  as  shown  in  Fig.  9*  Three  curves  are 
shown.  The  labels  1=1,6,19  indicate  the  number  of  i  terms 
in  the  expansion.  Normalized  a.,  values  determined  from  the 
derived  exact  formula  can  be  compared  in  Table  1  Wjjth  the 
values  of  a. .  determined  by  using  the  Hoff,  et  air  approxi¬ 
mate  formula1^  - 10) .  For  large  values  of  the  circumferential 
mode  number  n  these  normalized  values  are  nearly  constant  for 
all  values  of  n,  and  thus  essentially  independent  of  n.  It 
can  be  seen  that  the  values  obtained  from  the  approximate 
formula  are  very  good  approximations  to  the  exact  values  for 
n  greater  than  7  or  8.  At  the  lower  values  of  n,  there  is 
an  appreciable  difference. 


3.  Honeycomb  Cell  Shape  Equations 


As  indicated  previously  the  buckle  shape  of  crushed 
honeycomb  cells  is  to  be  approximated  by  the  modified  n=3 
radial  deflection  shape.  A  Fourier  series  which  will  repre 
sent  the  shape  is 
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where  b  is  half  the  effective  glue  line  width.  The  expression 
for  the  coefficient  a^^  is 
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Exact  I  Values 
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Pig.  9  Radial  Displacement  at  6«0  for  Different  Values  of  I 
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Table  1  Comparison  of  Fourier  Coefficients 


The  values  of  a.,.,  coefficients  have  been  computed.  To 
test  their  validity  tne^deflection  shape  produced  by  using  the 
a.  coefficients  for  the  degenerate  case  of  zero  width  glue 
line  has  been  compared  to  the  displacement  shape  produced  by 
the  a. .  coefficients  in  the  n=3  mode.  The  two  displacements 
agreed 

The  net  displacement  of  a  buckled  honeycomb  cell  is 
equal  to  the  Pourier  series  expression  for  the  displacements 
for  the  modified  n=3  buckling  shape  minus  the  original  hexa¬ 
gonal  cross  sectional  shane  which  is  the  n=6,  m=0  mode,  r.ince 
the  original  hexagonal  shape  is  in  the  m=0  mode,  the  displace¬ 
ment  shape  is  not  a  function  of  x. 

The  original  hexagonal  cross  sectional  shane  cylinder 
has  non  zero  i  and  j  terms  when  expanded  in  6j  9 ,  but  when  the 
hexagonal  shape  is  expanded  in  3J 0  the  odd  .1  terms  enual  zero, 
as  shown  in  Section  4  of  the  Appendix.  The  net  deflection  Is 
expressed  in  the  form 
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4.  Buckling  Stress 

As  indicated  in  the  previous  section  the  buckling  stress 
is  determined  by  using  the  expression  for  the  net  deflection 
to  compute  bending  and  stretching  strain  energy.  A  computer 
program  was  written  to  determine  the  strain  energies  and  the 
external  work  done  by  the  applied  load.  These  values  were 
then  used  to  compute  the  minimum  buckling  stress,  as  shown  in 
Section  5  of  the  Appendix .Buckling  stress  increases  with  cell 
wall  thickness  and  increased  effective  glue  line  width  as 
shown  in  Pig.  10.  In  this  figure  the  buckling  stress  is 
shown  in  non-dimensional  form  by  dividing  by  Young's  modulus  E. 
Cell  wall  thickness  t  is  non-dimensionalized  by  dividing  by 
the  cell  wall  width  c-  Note  that  the  buckling  stress  increases 
very  rapidly  with  thickness  for  thin  cell  walls  (assume  c  is 
constant)  at  small  values  of  t/c,  and  much  less  rapidly  after 
t/c  exceeds  0.02.  The  analysis  which  has  been  provided  is 
most  applicable  for  thin  wall  cells  and  becomes  less  accurate 
as  the  wall  thickness  increases.  This  is  because  the  cells 
will  tend  to  buckle  in  modes  which  do  not  represent  developable 
surfaces  as  the  wall  thickness  increases.  The  difference  in 
the  buckling  patterns  is  readily  observed  by  loading  cylin¬ 
drical  shells  of  different  wall  thicknesses  until  they  buckle. 
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Note  also,  in  Fig.  10,  the  significant  increase  in  buckling 
stress  as  the  glue  line  width  is  increased.  The  same  effect 
is  produced  by  decreasing  the  cell  size.  The  results  in  Fig.  10 
are  for  one  axial  buckle  length  (  1/ c  =  1.0).  However, 
increasing  the  glue  line  width  increases  the  buckling  stress 
significantly  at  other  values  of  the  buckle  length  (£  /c). 

This  is  shown  in  Fig.  11.  The  upper  group  of  curves  in  Fig.  11 
indicate  that  for  a  constant  effective  glue  line  width,  the 
buckling  stress  increases  as  the  value  of  £/c  decreases,  or 
as  the  buckle  length  becomes  shorter  and  shorter.  However, 
the  lower  group  of  curves  (  i  /c  =  4.0  and  i  / c  -  10)  show  a 
reversal  of  this  trend.  Although  this  behavior  seems  at  first 
glance  contradictory  it  is  in  fact  what  one  might  exoect. 

There  is  a  buckle  length  at  which  the  buckling  stress  will  have 
a  minimum  value  and  this  is  the  buckle  length  which  will  be 
observed  no  matter  what  the  cell  length  is.  If  this  were 
not  the  case  the  buckle  length  for  minimum  buckling  scress 
would  always  be  determined  by  the  overall  length  of  the  cell. 

It  has  often  been  observed  that  if  the  loading  of  a  cell 
is  stopped  soon  enough  a  number  of  folds  (with  a  short  buckle 
length)  will  be  formed  somewhere  along  the  cell,  usually  near 
one  end  or  the  other,  and  the  rest  of  the  cell  will  be 
unbuckled.  This  phenomenon  has  been  observed  for  both  static 
and  dynamic  loading.  The  results  in  Fig.  11  were  conrouted  by 
assuming  that  the  loaded  cell  would  buckle  at  some  arbitrary 
value  of  £/c.  Hence  these  results  should  not  be  interpreted 
to  mean  that  buckling  will  occur  at  any  of  these  values  of  1/ c 
shown.  It  can  be  said,  however,  that  for  t/c  =  0.001  if  an 
axial  load  is  applied  and  slowly  increased  until  buckling 
occurs  the  buckling  length  will  be  bounded  by  the  Inequality 
1  <  i/c  <10.  .  It  should  be  noted  that  for  80-0  -  1/2 

honeycomb  the  value  of  c  is  approximately  0.25in,  and  t  is 
approximately  0.007in.  Thus  t/c  =  0.028  and  for  £/c  =  4, 
jFl.Oin.  This  value  of  i  ,  the  buckle  length,  is  aDDreciably 
greater  than  the  values,  that  have  been  observed.  A  possible 
explanation  of  the  discrepancy  might  be  that  the  preferred 
buckle  length  (minimum  buckling  stress)  is  sensitive  to  the 
value  of  t/c,  and  the  results  would  therefore  have  been  con¬ 
siderably  different  had  t/c  =  0.028  been  used  in  the  compu¬ 
tations.  A  more  likely  explanation  Is  that  the  discrepancy 
is  simply  a  result  of  extrapolating  from  an  idealized  situa¬ 
tion  in  which  the  honeycomb  cell  Is  uniform  in  all  respects 
and  is  assumed  to  behave  elastically  right  up  to  the  time  of 
buckling,  to  the  real  situation  in  which  the  cell  is  far  from 
uniform  in  all  respects  and  buckling  is  more  plastic  than 
elastic.  The  results  should  be  interpreted  as  simply  an 
indication  of  how  the  buckling  stress  is  affected  by  the 
various  parameters,  rather  than  an  indication  of  the  buckling 
stresses  and  buckling  lengths  that  might  be  expected  under  a 
given  set  of  circumstances. 
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As  indicated  previously  a  constrained  axially  loaded 
hexagonal  cell  in  a  honeycomb  pad  buckles  in  a  modified  form 
of  the  n=3  mode  shape.  There  are  six  buckles  in  the  circum¬ 
ferential  direction.  Except  at  the  axial  nodes,  three  of  the 
buckled  sides  each  contain  two  of  the  straightened  out  axial 
folds  produced  in  the  construction  of  the  original  hexagonal 
cell.  The  maximum  number  of  buckles  in  the  axial  direction 
can  be  estimated  from  the  L/a  ratio  using  the  following 
expression: 


(A+Q)  <  , 
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ST) 
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(A+Q) 


m  <  /.  68  ~ 
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This  expression  is  obtained  by  considering  the  1! rat¬ 
ing  value  of  the  sine  of  the  angle  formed  by  the  buckled 
cell  wall  panel  and  the  original  cell  wall. 

For  comparison  it  nav  be  noted  that  the  crushed  hcr.c--- 
comb  cell  in  Fig.  6  has  an 


Hence  T 

m  £  1.68  -  =  21 

a 


From  the  photograph  about  19  buckles  can  be  counted  ale nr  the 
cell  length. 

Poisson's  ratio  for  paper  is  not  very  well  known.  To 
determine  it s,  influence  on  the  buckling  stress  of  paper 
honeycomb  the  computational  results  shown  in  Fig.  12  were 
obtained.  Obviously  Poisson's  ratio  has  an  insignificant 
effect  on  the  buckling  stress. 
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1.  Hexagonal  paper  honeycomb  cells  buckle  under  axial 
loading  into  a  pattern  which  approaches  a  developable  sur¬ 
face  for  very  thin  walls,  and  departs  more  and  more  from  that 
type  of  surface  as  the  thickness  of  the  cell  walls  increase. 

2.  The  effective  glue  line  width  between  adjacent 
honeycomb  cells  has  a  significant  effect  on  the  buckling 
stress  of  thin  wall  honeycomb  cells. 

3-  Variations  in  the  effective  glue  line  width  have 
a  greater  effect  on  the  buckling  stress  when  the  effective 
width  is  near  the  maximum  value  it  can  have,  i.e.,  the  cell 
wall  width. 

4.  The  cell  buckling  stress  increases  very  rapidly 
with  wall  thickness  at  thicknesses  less  than  0.02c  and  much 
less  rapidly  at  greater  wall  thicknesses. 

5.  The  buckling  stress  increase  with  wall  thickness 
becomes  more  pronounced  as  the  effective  glue  line  width 
increases . 

6.  Changes  in  Poisson's  ratio  have  an  insignificant 
effect  on  thQ  buckling  stress. 

7.  Honeycomb  pads  containing  hexagonal  cells  constrained 
by  other  hexagonal  cells  buckle  in  a  modified  form  of  the 
circumferential  mode  number  n=3  of  a  circular  cylindrical 
shell. 

8.  The  approximate  Fourier  coefficients  in  the  exoress- 
ions  derived  for  the  radial  deflections  by  Hoff,  Madsen  and 
Mayers  are  quite  good  for  surfaces  with  a  large  number  of 
circumferential  buckles  but  rather  poor  when  the  number  of 
circumferential  buckles  is  small. 

9.  A  more  precise  analysis  of  the  buckling  of  hexa¬ 
gonal  cells  is  not  justifiable  in  view  of  the  great  differ¬ 
ences  between  the  characteristics  of  an  actual  cell  and  the 
theoretical  cell  used  in  analysis. 


Recommendations 

To  obtain  a  paper  honeycomb  which  has  uniform  crushing 
strength,  special  attention  should  be  given  to  quality 
control  of  the  following  parameters: 

1.  Glue  line  width 

2.  Cell  wall  widths 

3.  Paper  thickness 

To  provide  further  authentication  for  this  analysis 
computations  should  be  made  of  buckling  stresses,  for  addi¬ 
tional  values  of  a/c  and  t/c.  Further  carefully  controlled 
experimentation  should  be  done  using  paper  honeycomb  with  a 
wider  range  of  values  of  c  and  t. 
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LIST  OP  SYMBOLS 
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a 


a 

A 


C^x) 

c2(x) 

E 

F 

h 

i 

j 

k 

% 

l 

x 

£y 

L 

m 

n 

P 

Q 

r 

t 

T 

u,  rr 


radius . 

(x)  Fourier  coefficient  of  .Fourier  series ' expansion  in  0. 

Fourier  coefficient  of  double  Fourier  series  expansion. 

Fourier  coefficient  for  n=3  mbde. 

aii6  Four-^er  coefficients,  for  n=6  mode  expanded  'in  6 jo  and 
A  3j9  respectively. 

radius  vector. 

maximum  inward  radial  di splacemeni:,  also  a  dummy 
variable .  * 

!  2  1 

Fourier  coefficient  normalized  with  respect  to  n  /a. 
half  width  of  effective  glue  line.  . 
dummy  variable.  . 

*  i  s 

buckled  honeycomb  cell  axial  half  :wave  length,  also 
cell  wall, width. 

radial  displacement  at  0  =  1 

radial  displacement  at  0  =  0.  i 

variable  equal  to  1  +  jn,  YoUng’s  modulus, 
variable  equal  to  1  -  jn.' 
thickness . 

axial  mode  number.  •  , 

circumferential' mode  dumber, 
variable  subscript. 

.  t  •  I 

length  > 

,  ,  l 

buckled  panel  axial  length. 

one  half  buckled  panel  circumferential  length, 
cylinder  length.  (  ! 

axial  buckle  mode  number. 

circumferential  buckle  mode  number. 

;  ! 

axial  load,  dummy  variable,. 

maximum  outward  radial  displacement. 

1 

radius,  radial  coordinate, 
thi'ckness.;  ‘ 

dummy  variable. 

axial  displacement,  axial  displacement  vector. 
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>1  In 


V  ue' 

Utot 

V,  V 


w,  w 


w 


W, 


net 

H 


n=3 

w 

n=3  mode 

x 

xk 

X 


Y 

6 

6m 

m 

ex' 


0  ' 


0 

0 


v 

e 

p 


* 


bending,  extensional  and  total  strain  energies. 

tangential  displacement,  tangential  displacement 
vector. 

radial  displacement,  displacement  vector. 

net  radial  displacement  of  buckled  hexagonal  cell. 

radial  displacement  required  to  buckle  circular 
cylinder  into  hexagonal  cylinder. 

radial  displacement  in  the  n=3  mode. 

radial  displacement  of  honeycomb  cell  in  modified 
n=3  mode. 

axial  coordinate. 

axial  coordinate  of  a  particular  cross  section, 
general  axial  location  of  axial  cross  sections, 
tangential  coordinate. 

position  vector  in  final  coordinate  system. 

angle  locating  slope  discontinuity  in  final  coordinate 
system. 

circumferential  period  =  2ir/n. 

total  axial  displacement  of  buckled  cylinders. 

Kronecker  delta. 

axial,  tangential  and  shear  strains  respectively. 

initial  position  polar  coordinate  at  x  =  X. 

initial  position  polar  coordinate  =  n/n  -  , 

at  x  =  -  X. 

angles  locating  slope  discontinuity  in  initial 
coordinate  system©  and  01#  respectively. 

Poissons  ratio. 

final  position  coordinate  system, 
mass  density. 

buckling  stress  of  honeycomb  cell. 

angle  between  axial  generator  and  buckled  panel. 

denotes  differentiation  with  respect  to  x. 


a 
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•APPENDIX 


DERIVATION  OP  EQUATIONS 


The  analysis  of  the  bucking  of  a  hexagonal  honeycomb 
cell  during  axial  loading  requires  five  consecutive  steps. 
These  are: 

1.  Determination  of  the  exact  displacements  which 
transform  a  circular  cylinder  into  a  polyhedral  developable 
surface. 

2.  Deriving  a  formula  for  the  coefficients  of  a 
Fourier  series  for  the  radial  deflection  expanded  in  terms 
of  the  circumferential  coordinate  0. 

3.  Derivation  of  a  formula  for  the  coefficients  of  a 
double  Fourier  series  for  expressing  the  radial  deflection 
in  terms  of  the  circumferential  coordinate  0  and  axial 
coordinate  x. 

4.  Derivation  of  a  formula  for  the  coefficients  of  a 
double  Fourier  series  for  including  in  the  radial  deflection 
of  buckled  honeycomb  cells  the  effect  of  the  effective  glue 
line  width. 

5.  Computation  of  the  strain  energy  in  cells  buckled 

in  a  specified  mode  and  then  using  this  energy  to  compute  the 
buckling  stress  for  that  mode. 

These  steps  will  now  be  discussed  in  detail. 

1.  Exact  Displacements 

The  initial  coordinate  system  of  a  circular  cylinder 
is  shown  in  Fig.  3.  A  circular  cylindrical  surface  can 
buckle- into  a  develooable  Dolyhedral  surface  such  as  the  one 
shown  in  Fig.  3  which  has  a  circumferential  mode  number  n=6 
and  longitudinal  mode  number  m*8.  In  this  polyhedral  devel¬ 
opable  shape  there  is  no  extensional  strain  at  any  point  in 
the  surface. 

In  order  to  determine  the  exact  radial  and  tangential 
displacements  for  a  circular  cylinder  that  has  buckled  into 
a  developable  surface  it  is  necessary  to  determine  the 
radial  displacement  amplitudes  as  a  function  of  x. 

Fig.  4  shows  orofile  views  along  the  longitudinal  axis 
of  a  general  developable  surface  at  angular  coordinate  loca¬ 
tions  0  =  0  and  0  =  ir/n.  The  radial  deflection  amplitudes 
as  a  function  of  x  are  determined  in  the  following  manner: 

Consider  the  longitudinal  cross  section  shape  at  0  =  0 
shown  in  Fig.  4a.  Coppa,  from  geometrical  considerations 
and  the  condition  of  no  extensional  strain,  has  shown  that 
the  radial  deflection  maxima  are: 


Preceding  page 
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Sin  21 
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(outward  deflection) 


/l  = 


a. 

Sin  22 


f**{ -%<***} 


(1.1) 


(inward  deflection) 


These  expressions  apply  only  if  the  cross  section  is  an  eoual 
sided  polygon,  which  it  will  be  at  the  cross  sections  where 
maximum  deflections  occur.  For  intermediate  sections  let  c. 
be  the  radial  displacement  at  0  =  0. 

Consider  0  <  x  <  i  . 

From  Fig.  4a,  c2(0)  =  -A,  c2(  i^)  =  Q. 

,  ,  _  ,  A  +  Q  \ 

cz  (x)  *  C— ■ — - )x-A  (1  2) 

From  Fig.  4b,  c^(0)  =  Q,  c^(  A  x)  =  -A. 

C,(x)s-  (  x  +  Q  -  Q-A  -  Ce  (x)  -  Cz  (j(x  -  x)  (1.3) 


Also  Cz(x)~  C,  (Jx~x) 


(1.4) 


Furthermore,  from  Fig.  4  it  can  be  seen  that  c2(x)  and  c,(x) 
are  periodic  in  x  with  a  period  equal  to  2iy  .  Fig.  A1  shows 
some  cross  sections  of  the  developable  surface  between  the 
axial  nodes.  Detailed  views  of  the  cross  sections  are  shown 
in  Figs .2a  and  2b.  Note  that  the  cross  sections  have  2n 
variable  length  sides  at  any  location  between  the  axial  nodes 

for  m  =  1,2, . ,M.  There  are  six  sides  corresponding  to  the 

circumferential  mode  number  only  at  the  axial  nodes,  x  =  m  i  , 
where  six  of  the  sides  degenerate  to  zero  length,  and  the  six 
remaining  sides  are  equal  to  2  l  or  2-n— ,  since  there  is  no 
extension  from  the  original  surface.  Axso  note  that  the  sides 
degenerate  to  zero  length  in  such  a  manner  that  it  aopears  that 
the  cross  section  at  an  axial  node  is  the  same  as  an  adjacent 
node  cross  section,  except  that  the  sections  are  rotated 
Vn  radians  with  respect  to  each  other. 
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Consider  now  a  general  cross  section  of  a  cylinder  buck¬ 
ling  in  an  arbitrary  circumferential  mode  n  as  shown  in  Pig.  A3. 
This  is  an  exaggerated  view  of  a  portion  of  a  cross  section 
between  the  axial  nodes,  such  as  the  cross  section  shown  in 
the  middle  of  Pig.  A2b. 

The  radial  and  tangential  displacements,  w  and  v,  shown 


change  the  original  cylindrical  surface  into  a  polyhedral 
developable  shape.  In  the  initial  unbuckled  coordinate  system 
0  is  the  angular  coordinate, and  £  is  the  angular  coordinate 
in  the  final  buckled  configuration  coordinate  system.  The 
angular  coordinate  locating  the  limiting  edge  of  one  of  the 
sides  which  is  a  function  of  x,  is  given  by  the  angle  8(x), 
measured  in  the  final  buckled  coordinate  system. 

Consider  first  the  deflection  for  £<B(x)  for  any 
general  position  x  =  X  along  the  cylinder  axis  as  shown  in 
Fig.  A3b. 

The  initial  position  point  0  =  e  ,  *  from  which  the 
final  position  point  £  =  8(x)  is  mappea'after  deflection  can 
be  determined  as  follows: 


at  £  =  8(x), 


e  =  eB(x) 


and  tan  £  =  tan  8 


For  inextensibility 


tan  {3  = 


a  +  c2  Cx) 


0  =  (f  1  c£ 

SL 


(1.5) 


From  Fig.  A3b 


tan  p  =  - 


3('-fx ) 


( /-  O  - 


+  Ce  f  x)  -  -for  f  £  X<  2J/X  ) 


(1.6) 


40 


Tangential  Displacement  s 


Consider  0  £  x  -  & 

A 

Substituting  Eq.  (1.6)in  Eq.  (1.5) 


_  [  &  +c2fx)J  *y  ( 

~ '  -  r  -  . 


e,M 

si  ( 0)  +  (6)  •+■ 


[a  +  cz  (*)  ] 
V(9)  =  z(f) 


*lLL£} 


Consider  the  deflections  for  5  ^  fi,  that  is, 


for 


J^vert  projection  =  a  cos  9  +  v(0)sln  9  +  w (&)cos6 

*  z  (f)  cos  f 

J^hor  projection  =  3.  sin  9  -  \/( ^cos  O  +  yv( S^sin  O 

=  [z  +  Cf  (x)  ]  tan  f 
Also  for  inextensibility 

[a.  •+  c2  (x)  2  t&n  ^  =  a.  9 

Substituting  Eq.  (1.9)  into  Eq.  (1.8) 


R  sin  O'  v(O)  cos  O  +  !a/  (9)  sin  O  =  el  9 

Multiplying  Eq.  (1.7)  by  cos  0  and  Eq.  (1.10)  by  sin  0  and 
adding  yields 


w  (9)  =  [a  +c2  ( x )]  cos  9  ■+  a.  6  sin  G  -  a. 


From  Eq.  (1.10)  and  (1.11) 

v  (9)  =  [a  ■+  c2  (x)]  sin  0  -  a  9  cos  6 
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(1.7) 


,  J.  •  O  / 


(1.9) 


(1.10) 


(1.11) 


(1.12) 


Note  that 


v(Q) 

Also  note 


-  +  a,  s/n  O 

de> 


(1.13) 


w(e)  -  +2LCos  9 - a.  -  2  a.  s/n  2  $ 

uf?  «“ 


Continuity  in  the  circumferential  direction  at  the 
deformed  cylinder  ends  demands  that ' *  : 


d  v 


dy  •- 


~ZTT 

ir  de  =  o 
Z9 


(1.14) 


Eq.  1.12  satisfies  this  requirement. 

The  deflections  were  determined  from  Pig.  A3a  in  a 
similar  manner  for  B(x).£5£l-  and  were  found  to  be: 

v/(x, e)  =  [a  *c,(x)]  cos(yj-6)*-(%-8)»sin(  % ~9)-a  cos(Z-6) 


M[,-fpe<-3,  U.: 


-  n 
^ ; 


v(x,e)=- 


=  -[a+c,  (x)] s/n  (  &)+(% -0)  a  cos  C%-&J 


as//?  C %-e) 

<t  6 


(1.16) 


Eqs .  (1.11),  (1.12),  (1.15),  and  (1.16)  together  comorise 
the  exact  radial  and  tangential  displacements  required  for 
a  circular  cylindrical  surface  to  form  into  a  developable  poly¬ 
hedral  shaped  surface.  The  terminology  "exact”  is  used  to  dis¬ 
tinguish  these  equations  from  the  series  approximation  forms 
which  are  subsequently  used. 


^3 


The  displacement,  u,  of  the  developable  surface  in  the 
axial  direction  is  determined  in  the  following  manner.  Refer¬ 
ring  to  Pig.  A4, 


u  -  -  x  ( /-  cos  0)  =  -x  (/  ~  - s/n 2 0  )  (1* 


where 

S/n  <p  =  ( A  +  Q) 

Prom  Eq.  (1.1) 

s/n  p  =  tan  ~ 


Inserting  Eq.  (1.19)  in  Eq.  (1.18)  yields 


2.  Fourier  Series  Expansion  in  0 

The  radial  deflection  given  by  Ecj.  (1.11)  can  be  exoanded 
into  a  Fourier  Series  in  terms,  .of  8  in  the  following  manner. 
Since  the  radial  deflection  is  periodic  in  0  with  period  2ir/n 
only  the  radial  deflection  between  the  circumferential  nodes 
is  required  for  the  Fourier  Series  expansion.  The  radial 
deflection  at  any  axial  location  X  can  be  expressed  as: 


oo 


*/  ( X,  0)  =  £  3-j  ( x)  cosjn 9+  loft) s/n jn 6  +  aD(x) 
jn 


(2.1) 


and 


Since  the  radial  deflection  is  an  even  function  in  e,  b  =  0 

J 


r 

2 


C /_  ?J) "(x>  e)cosjne‘,e 


(2.2) 


:x 


From  symmetry 


(2. 


ij(x>  =rf  (/~i  O  w  (y> e)  c°sjnede 


e,(x)  o  g 

-  Cfo  (f~i  fj  )co$jr>0d9*J'  ( '-LSj  )w(X,e>d^Cosjn9J^{  2 . 


where 


*,(*)-?  U-frJ  -  %[,-*£] 


(2. 


Now  considering  the  symmetry  of  the  deformation  (Pig.  H) 


w (X)  -for  e(X)>  Qfi  (x)  -  W(JA~X)  for  0,  f/x~ X)  *  6^,  fc/Ji 2.' 


And  since  it  is  also  symmetrical  about  x 
must  be  considered. 

When 


y  then  Of(Jx-X) 


*x,  only  0  X 

'  =  (A-x) 

=  O 


l 


x 


(2.1 


where  0,(4  -X) 

x  =  £  -X 
x 


is  the  initial  rosition  angular  coordinate  at 


*/  d -*)m  £  L'~ 


Jyf  /777TX 

a/x  t}L  (2.6 


\nd  0qi(ax-X)  is  the  slope  discontinuity  angular  linut  in  the 
initial  position  angular  coordinate  at  x  =  z  -  X. 

To  facilitate  reduction  of  the  second  integral  in  Eo.  (2  H) 
transform  it  to  the  variable  e1(4x-X). 


.*/  (J*-*)  “i-o(x) 


I 


9(x)  =  X-e,  (J, - x)  =  %-e,  (4-1O 


(2.10) 


de(x)  =  -  do,  Ux-X) 

Substituting  Eqs.  (2.6)  through  (2.l\3)  into  Eq .  (2.t)  yields 

*  re*M  s  . 

.aj00-7?lj  ) w(X ,0)  cosjna/e 

+  0  f  _  (/-^j)^C 4  -X,  3‘J  CO  5//?  r /-  Y'dQ,)] 

?  A r0g(x) 

slj  (X)-  ~  ('-$ tj)  [  j^(X,e)  cosjn  c/e 

r&  uasf  *o  1 

+J  *  £(j, -X,  J  cos j n  (Z-9\)dQ,]  (2.11) 


C 05 jn  (Y~0,)  =  COS jn(%~ 9,)  =  COS  (j*  yo9,)=(~ l)Jcosjr  6,'  (2.12) 


and  from  Eq..  (1.11) 


"(Jx-X,  9,) =  [i  +  C$(j?x-X)]cdsG,  '+  &  9 /Sin  Of -a. 


but  from  Eq.  .(1.3) 


cl(Jx-x)~  c,(X)  o<X<!> 


(2.13) 


(2. 1^) 


Evaluating  the  first  integral  in  Eq.  (2.11) ! by  use  of  Eqs.  (1.11) 
yields  < 

,  .  Cr  '  x->  /  r  Sin £9  sinF01  ^rsioEB 

J  *v (X,  9) cosjn9c/9 s  f fa  +  Cg  (X)j g [  ~JT~  +  £& 

O  '•  *v 

0*0 

,  sin  F9 7  30r  cosEe '  COsFOn)  ,  a.S/njn80 

=  “T2“ 2  (-  ~T~  F  JJ  ~  Jn 

•  1 

hi 

F  =  /  +  )n  F-  /-in 


t7 


. 


Eq.  (2.17)  reduces  to 


V 


+  Sir*  %( 


cosEO^i  cosF6fl 


) 


3.  Determination  of  Fourier  Coefficients  for  Double  Fourier 

Series 

The  expression  given  in  Eq.  (2.2)  used  with  the  cerived 
Fourier  coefficient  in  Eq.  (2.18)  yields  a  deflection  st.ape 
which  has  no  slope  discontinuities  in  the  circumferential 
direction  but  has  slooe  discontinuities  in  the  axial  direction. 
This  is  because  a.(x)  has  slope  discontinuities  in  the  axial 
direction.  In  order  to  obtain  the  radial  displacement  func¬ 
tion  which  has  no  axial  direction  slope  discontinuity  for  a 
finite  number  of  terms,  the  radial  deflection  can  be  expressed 
as  a  double  Fourier  series  which  is  the  same  as  expanding 
a^(x)  in  a  Fourier  series  as  a  function  of  x  or: 


OO  CO 

W(X,9)  -£  Y,  */j  C0S  C0SJn 

i-o  j~0 


C 


Using  aj(x)  as  it  is  defined  in  Eq.  (2.1) 

m 

*ij  -  2~ f  )  *j(*)  COS  '-ZLH*  dx 

-  L 


m 

For  a^  an  even  function  in  x 

L 


*/j  =  *?/  [0-i*i)*j Mje»s  lal*  Jx 


(3.1) 


(3.2) 


(3-3) 
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When  Eq.  (2.18)  is  substituted 


Substitution  of  Eqs.  (3.5)  and  (3.6)  into  (3.4)  after  simplification  yields: 


h .  Fourier  Coefficients  for  Honeycomb  Cell 

Inspection  of  the  buckled  cells  of  the  impacted  paper 
honeycomb  specimens.  Fig.  6,  shows  that  the  radial  displace¬ 
ments  from  the  sides  which  form  the  original  hexagonal  shape 
are  alternately  inward  and  outward  as  one  progresses  around 
the  cell  wall.  The  original  hexagonal  shape  can  be  represented 
by  the  analytical  expression  for  the  radial  displacement  from 
the  derived  Eqs.  (3.8)  and  (3.1)  by  using  the  longitudinal  mode 
number  m  =  0  and  circumferential  mode  number  n  =  6. 

Cross  sections  of  a  buckled  cylindrical  shell  in  the 
n  =  3  mode  have  three  sides  which  form  an  equilateral  tri¬ 
angle  at  the  axial  nodes.  However,  in  the  n=3  mode,  there  are 
6  sides  (2n)  at  any  cross  section  between  the  axial  nodes. 
Maximum  displacements  in  the  buckled  honeycomb  occur  at  the 
circumferential  mid-points  of  the  original  sides  of  the 
hexagon.  The  difference  between  the  maximum  inward  deflection 
at  these  points  for  the  n  =  3  mode  and  the  n  =  6,  m  =  0  mode 
is  equal  to  the  maximum  outward  deflection.  Fig.  7.  Thus 
the  alternating  inward  and  outward  displacements  observed  in 
Fig.  6  are  the  same  as  the  displacements  of  a  circular  cylin¬ 
drical  surface  in  a  modified  n  =  3  buckling  mode  minus  the 
n  =  6,  m  =  0  mode  which  represents  the  original  hexagonal 
cylindrical  shape  of  the  cell. 

The  modified  n  =  3  buckling  mode  of  the  honeycomb  cells 
is  a  function  of  the  width  of  the  glue  line  which  remains 
intact  during  the  buckling,  i.e.  "the  effective  glue  line 
width,'*  see  Fig.  A5.  Here  b  is  the  half  width  of  the  effec¬ 
tive  glue  line.  The  cross  section  of  a  buckled  circular 
cylindrical  shell  at  the  midpoint  between  the  axial  nodes  of 
the  n  =  3  buckling  mode  is  the  same  as  the  cross  section  for 
the  n  -  6,  m  =  0  mode.  This  means  that  the  deflections  are 
the  same  and  consequently  there  is  no  net  displacement  at 
that  section  between  the  two  deformation  modes.  This  axial 
location  corresponds  to  the  end  conditions  of  a  honeycomb 
pad  in  Fig.  A6  cell,  where  displacements  are  constrained  by 
the  facing  paper.  The  radial  displacements  of  the  buckled 
honeycomb  are  determined  in  the  following  manner: 
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Pig.  A5  Modified  Form  of  n  =  3  Mode  Shape 


expressed  a|;FlB'  ®’  the  radlal  deflecti°n  in  the  modified  n  =  3  shape  may. he 


Integrals  are  of  the  form: 
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5.  Strain  Energy  Expressions 

To  obtain  the  critical  buckling  stress  the  strain 
energy  is  calculated  and  then  equated  to  the  work  done  by  the 
force  required  to  cause  the  buckling.  The  strain  energy  cal¬ 
culation  proceeds  as  follows.  Transform  the  axial  coordinate 
x,  in  Eq.  (4.1)  to  a  new  axial  coordinate  x~  (see  Fig.  8). 
This  is  done  in  order  to  satisfy  the  natural  boundary  condi¬ 
tions  of  zero  deflection  and  simple  support  at  the  limits  of 
integration  in  the  strain  energy  expression.  The  transfor¬ 
mation  relationships  are: 

v  =  0  -  X*  -  O  at  the  limits 

y  ss  v  -  +  ..A.*  between  the  limits 

Z  TTX 


Insert  this  expression  in  Eqs.  (3*1)  and  (4.1)  and 
subtract  to  obtain  the  net  radial  deflection  of  the  hexaro^i 
cell.  Note  that  w(x5e)  when  n  =  6,  m  =  0  reduces  to 

eo  oo 

£  £  a-ijt  cos  6/0 

i-O  j-0 

which  can  be  written 


/  * 
~1 


Cos  3jO 
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This  reduces  Eq.  (5.5)  and  (5.6)  to 


third  and  higher  powers  of  derivatives  except  in  u.  produces 


u(x)  is  axial  displacement  function 
u(o)  =  0,  the  total  displacement  =  u(L) 

For  one  buckling  length  u  is  defined  in  Eo.  (1.20). 


“  (t’ih  ('-*&)) 


(5.15) 


(5.16) 


From  Eq.  (1.18)  it  can  be  seen  that  u  is  a  function  of  sin  0 
which  can  be  determined  in  the  following  mannev': 


6  7 


3  deflection  shape,  w  is  minimum  at  x  =  0,  and  maximum 


